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Abstract. In the young stellar object (YSO) phase of their lives, massive 
stars drive bi-polar molecular outflows. These outflows produce beautiful, often 
hourglass shaped, cavities. The central star possesses a powerful stellar wind 
(v ~ 2000 kms^ 1 ), and possibly a dense equatorial disk wind (v ~ 400 kms -1 ), 
which collide with the inner surface of the bi-polar cavity and produces hot 
(T ~ I0 5 — 10 8 K) shocked plasma. A reverse shock is formed at the point where 
the ram pressure between the preshock flow balances the thermal pressure of the 
postshock flow and provides a site for the acceleration of non-thermal particles 
to relativistic energies. Hydrodynamical models of the wind interaction, coupled 
with calculations of the non-thermal energy spectrum, are used to explore the 
observable synchrotron and gamma-ray emission from these objects. 



Introduction 



The formation of massive stars involve s outflows ( Garay Sz Lizanolll999l : Reipurth &z Bally 
120011 : iBaneriee k Pudritdl2006l . 120071') . Bi-polar cavities around YSOs are com- 
monly observed in star formation (|Garav & Lizandll999T). with collimation fac- 
tors of ~ 2 - 10 for massive YSOs fMYSOs) (|Beuther et alJ l2002fc iDavis et al. 
2004). Numerical simulations of high-mass star formation show that outflows 
are a by-product of the collapse process which for ms the massive star, and 
may be due to radiation p ressure and mag netic fields (|Yorke Sonnhaiteri r2002; 
iBanerjee fc Pudritd 2007). Th e termination shock of an MYSO jet can provide 
a site for particle acceleration ( Araudo et al .1 [20071 ) . 

The widths of IR recombination line emission observed from MYSOs in- 
dicates the presence of dense outflows with ve locities ranging from 100 to > 
340 kms -1 ( Drew et al. 1993 ; Bunn et al. 19951 ). Further confirmation of out- 



flows has come fro m high angular resolution radio observations (e.g. iHoare et al 



1994 ; Hoare 2006). One explanation would be a disk wind generated when UV 



flux from the star is absorbed and re-emitted by the material at the surface of 
the disk (e.g. iDrew et al.lll998T l. 



X-rays have been detected from deeply embedded MYSOs in star formin 



regions by the Chandra X-ray observatory (hereafter Chandra) (e.g. IWang et all 
20071 . 120091 ). An initial observation of Mon R2 detected X -ray emission from 



the intermediate mass objects consistent with the hard and highly time variable 
X-ray emission cau sed by magnetic flaring activity between the star and disk 
( Kohno et al. |2002| ). Further analysis of the same data, coupled with high res- 
olutio n near-IR interferome try identified the separate constituents of Mon R2 
IRS3 (jPreibisch et alJ 120021 ). The X-ray emission from IRS3 A and C, with a 
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count rate of 0.166 ± 0.041 ks 1 for the former, could not be explained by the 
standard scenario for massive stars (i.e. wind embe dded shocks produc ed by in- 



stabilities inherent in radiatively-driven winds - see Owocki et al.|[l988l ). yet the 
estimated stellar masses of these objects implies they will have radiative outer 
envelopes which poses problems for the generation of X-rays through magnetic 
recombination between the star and disk. In this work we examine another po- 
tential source of X-rays, namely the collision between the outflowing stellar and 
disk winds and the inf ailing envelope. 

We find that the interaction of the stellar and disk wind with the cavity wall 
produces shock heated plasma which emits X-rays in agreement with Chandra 
observations (Parkin et al. 2009, in prep). The reverse shock, plus weak shocks 
within the hot gas, provide a site for particle acceleration. Therefore, a timely 
question, considering the recent launch of the Fermi Gamma-ray Space Telescope 
(Fermi), is whether the winds-cavity interaction produces an observable non- 
thermal energy spectrum extending up to 7-ray energies. To this aim we have 
developed our existing model for thermal X-ray emission to consider non-thermal 
radio and 7-ray emission. 



2. Winds- cavity interaction 

The model consists of a MYSO situated at the centre of a previously evacuated 
bipolar cavity which is surrounded by infalling molecular cloud material. We 
include the stellar wind and a disk wind which emanates from the surface of the 
circumstellar accretion disk; both are assumed to be at terminal velocity. Due to 
the spatial scales under consideration we do not attempt to model the radiative 
driving of the win ds as this requires h i gh spatial resolut ion in the vicinity of 
the star /disk (e.g. lOwocki et al.l Il994l ; iProga et al.l Il998l ). For simplicity we 



adopt an angle dependant wind prescription, whereby the stellar wind occupies 
a region up to a polar angle of 60° and the disk wind occupies the region from 
60° — 85° . The density and velocity distibutions for the infalling cloud material 
are described by the equations of iTerebev et al.1 ( 1984! ). The morphology of 



the pre-existing outflow cavity i s desc ribed by a simple analytical prescription 
similar to that of I Alvarez et al.l ( 2004! ). 

We consider a ~ 30 M Q 08V star with a mass-loss rate of 10 -7 M Q yr _1 and 
terminal wind speed of 2000 kms -1 . The disk wind has a mass-loss rate of 
10~ 6 M Q yr _1 and a terminal wind speed of 400 kms -1 . The unshocked winds 
are assumed to be at a temperature of 10 4 K. The mass infall rate for the cloud 
is 2 x 10 _4 MQyr _1 and the cavity opening angle is 30°. The winds-cavity 
interaction was followed for a simulation time of t = 2000 yrs. 

Examining Fig. Q] we see that the disk wind, which has a post-shock temper- 
ature of ~ 10 6 K, lines the cavity wall beyond which resides the infalling cloud 
material. The hotter post-shock stellar wind fills a volume closer to the pole. 
A reverse shock is established within ~ 10 16 cm of the star/disk which can be 
clearly seen in the temperature plot of Fig. [TJ The shape and position of the re- 
verse shock fluctuates with time causing variability in the observed emission on 
timescales of ~ 50 yrs. The shocked stellar and disk winds reach temperatures 
up to ~ 10 8 K and ~ 10 6 K, respectively. 
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Figure 1. Density (left) and temperature (right) snapshots taken from a 
hydrodynamic model of the wind-cavity interaction at t = 1000 yrs. The tick 
marks on the axis correspond to a distance of 10 16 cm. 
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Figure 2. 1-100 GHz radio spectrum calculated from the hydrodynamic 
simulation shown in Fig. [T] using an inclination angle of 60°. In this plot 

CB = Crcl = 1%. 
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Figure 3. Same as Fig. [2] except £ B = 5. The contribution to the free-free 
flux from the unshocked stellar and disk winds and the shocked plasma are 
shown. 



Thermal X-ray emission in agreement with Chandra observatio ns of MYSOs 
(e.g. iKohno et afll2002t iPreibisch etaD 120021 : iGiardino et~al"1l2004l ) is produced 
by the shocked gas adjacent to the reverse shock. In particular, hard X-rays 
are produced (E > 5 keV), the simulated ACIS-I count rate is ~ 0.2 ks~\ and 
the visual extinction to the star is ~ 30 — 100 mag (corresponding to a column 
of 5.7 — 190 x 10 21 cm 2 ). We note that the emission weighted column density 
( Parkin &; Pittard 20081 ). which is indicative of the column density to the regions 
of highest intrinsic emission, is higher than the column density to the star. 



3. Radio Emission 

The reverse shock present in the simulations provides a site for the acceleration 
of non-thermal particles to relativisitic energies. We calculate the synchrotron 
emission and absorption, the Razin effect, and Inverse Compton (IC) cooling 
of electrons. The non-thermal and magnetic field energy densities are assumed 
to be proportional to the thermal energy density. The non-thermal particle 
distribution is described by a power-law distribution with p = 2. The ratio 
of the magnetic to thermal energy density is £b> and similarly the ratio of the 
relativistic particle to thermal energy density is Q e \. The free- free emission from 
the ionized winds is also calculated. 

Fig. [2] shows the resulting radio spectrum when (3 = Crei = 1% is assumed. 
The total spectrum has a flat slope, indicative of optically thin plasma. The 
free-free emission from the unshocked winds is lower than the contribution from 
the shocked winds inside the cavity. The thermal emission is dominant; the 
synchrotron emission only has a small impact on the total spectrum at low 
frequencies. Increasing the post-shock magnetic energy density (£3 = 5%) causes 
a rise in the magnitude of the synchrotron emission (Fig. [3]) and alters the slope 
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Figure 4. Non-thermal energy spectrum calculated from the hydrodynamic 
simulation shown in Fig. Q] where Cb = Crci-c = Cmi-p = 5%. The sensitivit y 
curve for the Chercnkov Telescope Array (CTA) is from Bcrnlohr ct al. (2008). 
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Figure 5. Same as Fig. |4] except Crci- P = 50% 
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of the oberved spectrum. With this in mind the observed radio spectrum could 
be used to constrain parameter values for the magnetic and relativistic particle 
energy density. However, a major difficulty with this approach lies in attaining 
observational constraints due to the low surface brightness and large spatial 
extent of the emission region; typical observations of massive star forming regions 
strive for high spatial resolution which resolves out emission on the scales we 
are interested in. 

4. Gamma-ray emission 

To examine the possibility of observable 7-ray emission we calculate the emission 
from: Relativistic Bremsstrahlung, IC emission, and Pion decay. Fig. [5] shows 
that IC emission is dominant at low energies whilst Pion decay is dominant at 
higher energies. For £b = Crei-e = Crei-p = 5% the observed flux is considerably 
weaker than the sensitivity of Fermi, HESS, HESS-II, and CTA. In the case of 
efficient particle acceleration (Crel-p = 50%) this gap is reduced and there is a 
tentative possibility of a detection with CTA (Fig. [5]). These results then imply 
that 7-ray emission will not be detectable from a single MYSO. 

We close with the note that the possability of detecting non-thermal 7-ray 
emission may be increased if consideration is given to the cumulative effects of 
a cluster of massive stars, or, in the case of efficient particle acceleration, from 
single MYSOs more massive than our test case 08V star. 
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